Nuclear collective excitations using correlated realistic interactions: 
the role of explicit RPA correlations 
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We examine to which extent correlated realistic nucleon-nucleon interactions, derived within the 
Unitary Correlation Operator Method (UCOM), can describe nuclear collective motion in the frame- 
work of first-order random-phase approximation (RPA). To this end we employ the correlated Ar- 
gonne V18 interaction in calculations within the so-called "Extended" RPA (ERPA) and investigate 
the response of closed-shell nuclei. The ERPA is a renormalized RPA version which considers explic- 
itly the depletion of the Fermi sea due to long-range correlations and thus allows us to examine how 
these affect the excitation spectra. It is found that the effect on the properties of giant resonances is 
rather small. Compared to the standard RPA, where excitations are built on top of the uncorrelated 
Hartree-Fock (IfF) ground state, their centroid energies decrease by up to 1 MeV, approximately, in 
the isovector channel. The isoscalar response is less affected in general. Thus, the disagreement be- 
tween our previous UCOM-based RPA calculations and the experimental data are to be attributed 
to other effects, mainly to a residual three-body force and higher-order configurations. Ground-state 
properties obtained within the ERPA are compared with corresponding HF and perturbation-theory 
results and are discussed as well. The ERPA formalism is presented in detail. 

PACS numbers: 24.30.Cz, 24.30.Gd, 13.75.Cs, 21.30.Fe, 21.60.-n, 21.60.Jz 
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I. INTRODUCTION 



For the description of nuclear collective excitations and 
especially giant resonances throughout the nuclear chart, 
the theory of choice is the random phase approximation 
(RPA) — and its various siblings: the quasi-particle RPA 
for open-shell nuclei, the renormalized RPA, the second- 
order RPA (SRPA) and of course the relativistic versions. 
The most primitive applications of the usual RPA — i.e., 
non-relativistic, first-order, for closed-shell nuclei — con- 
sist typically of a ground-state description based on the 
Woods-Saxon potential and a phenomenological resid- 
ual particle-hole (ph) interaction, e.g., of Landau-Migdal 
type. An important development in the field has been 
the construction of increasingly refined energy function- 
als and corresponding interactions, with those of Skyrme 
type dominating the landscape. These have allowed what 
is termed "self-consistent" RPA calculations: the residual 
ph interaction is derived from the same interaction used 
in the Hartree-Fock (HF) description of the ground state. 
Self-consistency ensures that certain important formal 
properties of the RPA solutions hold exactly, including 
the separation of spurious states and the preservation of 
sum rules. Moreover, physical links between the bulk 
properties of nuclear matter and those of collective nu- 
clear states can be established. The self-consistent RPA 
offers also enhanced predictive power: By starting from 
a more and more refined energy functional one wishes to 
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make reliable extrapolations towards the unknown ter- 
ritories of the nuclear chart. Phenomenology enters this 
scheme in the form of the effective two-nucleon (NN) and 
thrce-nucleon (NNN) interaction. Contrary to the bare 
NN and NNN interactions, these are tailored for the use 
in mean-field type calculations such as HF and RPA. An 
ansatz is typically made for the effective NN(N) force and 
an appropriate fitting strategy has to follow. 

Highly accurate parameterizations of the bare NN force 
are nowadays available At the same time, at- 

tempts are made to derive the NN, NNN, and many- 
nucleon force starting from first principles, within chiral 
perturbation theory The possibility is currently be- 

ing explored to combine mean-field theory with such re- 
alistic NN potentials. Two methods have been developed 
recently for "taming" the NN interaction and thus mak- 
ing it usable within models like HF and RPA: the Viow-fe 
low-momentum interaction is derived within renormal- 
ization group theory 6] and the Vucom correlated NN 
interaction is derived using the unitary correlation op- 
erator method (UCOM) lllllil- Although constructed 
following different formalisms, the two potentials have 
similar low-momentum matrix elements. In this work we 
focus on the UCOM potential. 

Within the UCOM, the major short-range correla- 
tions, induced by the strong repulsive core and the tensor 
part of the bare NN potential, are described by a state- 
independent unitary correlation operator. This can be 
used to introduce correlations into an uncorrelated many- 
body state or, alternatively, to perform a similarity trans- 
formation of an operator of interest. Applied to a realis- 
tic NN interaction, the method produces a "correlated" 
interaction, Vucom , which can be used as a universal ef- 
fective interaction, for calculations within simple Hilbert 
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spaces. 

The aim of the UCOM is to treat exphcitly only the 
state-independent short-range correlations; long-range 
correlations should be described by the model space. 
This tells us already that the UCOM-based HF is not 
enough, since a Slater-determinant wavefunction is un- 
able to describe correlations. It was found indeed that, 
although bound nuclei were obtained using the Vucom 
already at the HF level, the binding energies were un- 
derestimated by about 4 MeV per nucleon. Second-order 
perturbation theory, however, constitutes a tractable and 
adequate extension to the "zero-order" description pro- 
vided by HF, as far as nuclear binding energies are con- 
cerned p!o| . 

In a recent publication we employed the UCOM 
Hamiltonian in standard, self-consistent RPA calcula- 
tions to study nuclear giant resonances. The ground state 
was described by the uncorrected HF state, as usual. 
The main focus was on the isoscalar (IS) giant monopole 
resonance (GMR), the isovector (IV) giant dipole res- 
onance (GDR), and the IS giant quadrupole resonance 
(GQR). Highly collective states were indeed obtained for 
various closed-shell nuclei ranging from ^^O to ^"^Pb. We 
achieved a reasonable agreement with the experimental 
centroid energies of the IS GMR. By contrast, the ener- 
gies of the IV GDR and the IS GQR were overestimated 
by several MeV. 

Obviously, the Vucom is not a traditional effective in- 
teraction. Partly because no long-range correlations are 
(effectively) included in the UCOM, the corresponding 
nucleon effective mass in nuclear matter obtained in a 
HF calculation is very low (less than half the bare nu- 
cleon mass) . This is confirmed by the HF results in finite 
nuclei: the single-particle level density was found too low. 
It is also manifested by the above-mentioned RPA results 
on the GQR and GDR centroids. From the preceding 
discussion it follows that, besides the possible important 
role of missing three-body terms in the Hamiltonian (see 
Seem, another source of our failure to describe nuclear 
collective states quantitatively can be the inadequacy of 
the RPA method to take into account residual long-range 
correlations. 

The standard RPA is based on two assumptions which 
hint at possible remedies. First of all, only one-particle- 
one-hole excitations are taken into account. One can 
include higher-order configurations, starting with two- 
particle-two-hole within SRPA. Given that an extended 
model space is of great importance when using the 
Vucom, it is imperative to examine the effect; work along 
this line is in progress. The second assumption is that 
one can approximate the true RPA ground state by the 
HF ground state. One argument follows from the fact 
that the RPA is the theory of small-amplitude vibra- 
tions around the HF ground state, therefore correlations 
should be small anyway for the use of RPA to be justi- 
fied. It is not obvious that this assumption holds when 
the UCOM Hamiltonian is used, given the large correc- 
tion to the HF binding energies due to second-order 



and RPA correlations. Therefore, in this work we 
wish to examine the effect of explicit RPA ground-state 
correlations on our results for nuclear collective states. 

To this end we use a renormalized version of the RPA, 
developed in Refs. [l^ [14. .15], where excited states are 
built on top of the true RPA ground state. It is for- 
mulated in the single-particle basis which diagonalizes 
the one-body density matrix, i.e., the natural-orbital ba- 
sis, and its equations are solved iteratively. Following 
Ref. [Ill, we will call it Extended RPA (ERPA). The 
term has been used in the literature to describe also 
higher-order RPA models (such as SRPA) , but this is not 
the case here. Renormalized formulations of the RPA, 
such as the ERPA, consider explicitly the depletion of 
the Fermi sea in the ground state due to RPA correla- 
tions. Besides corrected excitation properties, the ERPA 
allows us to evaluate also corrected ground-state proper- 
ties, namely single-particle energies and occupation num- 
bers. It is derived using the number-operator method, 
which, contrary to the quasi-boson approximation, does 
not suffer from double-counting the second-order contri- 
butions [Tsl.!!^. 

This paper is organized as follows. In Sec.^we outline 
the UCOM method and introduce the correlated Hamil- 
tonian. In Sec. lllll we present the ERPA method and the 
formalism used to evaluate ground-state properties and 
transition strengths. In Sec. II VI we present our results. 
In Sec. we summarize and discuss future perspectives. 

II. THE CORRELATED HAMILTONIAN 

In this work we employ a correlated interaction con- 
structed within the unitary correlation operator method 
(UCOM), based on the Argonne V18 nucleon- nucleon 
(NN) interaction. Here we outline the basic principles 
of the UCOM scheme 0, S H [l3- More detailed yet 
compact descriptions of the method can be found in 
Refs. ESInl. 

The basic idea is the explicit treatment of the 
interaction-induced short-range central and tensor corre- 
lations. These are imprinted into an uncorrelated many- 
body state l^*) (e.g., a Slater determinant) through a 
state-independent unitary transformation defined by the 
unitary correlation operator C, resulting in a correlated 
state \^), 

1*) = ^!*). (1) 

The correlation operator C is written as a product of 
unitary operators Cq and Cr describing tensor and cen- 
tral correlations, respectively. Both are formulated as 
exponentials of a Hcrmitian generator, 

C = CnCr = exp[-i^5n_ij] exp[-i^5rr_y]. (2) 

The construction of the two-body generators gr and gji 
follows the physical mechanisms by which the interac- 
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tion induces central and tensor correlations. The short- 
range central correlations, caused by the repulsive core 
of the interaction, are introduced by a radial distance- 
dependent shift pushing nucleons apart from each other 
if they are within the range of the core. Tensor corre- 
lations between two nucleons are generated by a spatial 
shift perpendicular to the radial direction. For a given 
bare potential, the corresponding correlation functions 
are determined by an energy minimization in the two- 
body system for each {S, T) channel. 

Matrix elements of an operator O with correlated 
many-body states [vP) can be equivalently written as 
matrix elements of a "correlated" (transformed) opera- 
tor O = C^OC and uncorrelated many-body states |^), 
thanks to the unitarity of the correlation operator. Thus, 
one can work in simple Hilbert spaces (simple states) us- 
ing correlated operators, as well as with bare operators 
and explicitly correlated states. By applying the unitary 
transformation to a bare NN interaction, a phase-shift 
equivalent correlated interaction is obtained, which is 
suitable for use in tractable model spaces M,hOlhJ. The 
same transformation can then be applied to any other op- 
erator under study, as is needed for a consistent UCOM 
treatment. 

There are a couple of important issues arising in ac- 
tual applications. In an A-body system the correlated 
operator contains irreducible contributions to all particle 
numbers. Within a cluster expansion of the correlated 
operator 

O = C^OC = Ol^l + 0[21 + . . . + Ol^l , (3) 

where O^"] denotes the irreducible n-body contribution, 
we usually employ a two-body approximation, i.e., three- 
body and higher-order terms of the expansion are ne- 
glected. Starting from the uncorrelated Hamiltonian H 
for the A-body system, consisting of the kinetic energy 
operator T and a two-body potential V, the formalism 
of the UCOM is used to construct the correlated Hamil- 
tonian in two-body approximation 

ijC2 ^ +f[2l +t>[2l =T + VuCOM, (4) 

where the one-body contribution comes only from the 
uncorrelated kinetic energy T^-^^ — T. Two-body contri- 
butions arise from the correlated kinetic energy T^^] and 
the correlated potential V^'^\ which together constitute 
the phase-shift equivalent correlated interaction Vucom- 
It has been verified that higher-order contributions due 
to short-range central correlations can be neglected in 
the description of nuclear structure properties 0. The 
tensor interaction, on the other hand, is long-ranged and 
thus generates long-range correlations in an isolated two- 
nucleon system. However, the long-range tensor correla- 
tions between two nucleons embedded in a many-nucleon 
system are suppressed by the presence of other nucleons, 
leading to a screening of the tensor correlations at large 
interparticle distances. In terms of the cluster expan- 
sion this screening appears through significant higher- 
order contributions. In order to avoid large higher-order 



contributions, and at the same time to effectively de- 
scribe the screening effect, the range of the tensor correla- 
tion function — more precisely, the "correlation volume" 
^i^ '^'' 03 — is restricted during the parameterization 
procedure. Restricting the range of the tensor correla- 
tor has another important function, namely to ensure 
that only state-independent, short-range correlations are 
described by the UCOM. By varying the correlation vol- 
umes — the only parameters entering the formalism — 
a family of correlators and respective correlated interac- 
tions are obtained. 

The question is then how to optimize these parame- 
ters in order to best describe the screening effect and the 
separation of the two types of correlations. As demon- 
strated in Ref. 03 1 this can be done with the help of 
exact few-body calculations. In particular, applications 
within the no-core shell model showed that for the Ar- 
gonne V18 potential the value 4^'°^ = 0.09 fm^ leads 
to the best description of binding energies in ^H and 
''He. For this choice of tensor correlator range the miss- 
ing genuine three-nucleon interaction and the ommitted 
higher-order terms of the cluster expansion of the corre- 
lated Hamiltonian effectively cancel each other. As was 
subsequently shown within many-body perturbation the- 
ory and verified by RPA calculations 12], this can- 
celation remains at work throughout the nuclear chart, 
as far as the binding energy is concerned. The same does 
not necessarily hold, however, for other ground-state ob- 
servables, such as charge radii, or for excited-state prop- 
erties. 

In this work we will use the correlated Argonne V18 
potential with 4^'°^ — 0.09 fm'^. No tensor correlator 
is employed in the triplet-odd channel, where the tensor 
interaction is much weaker. We start from a Hamilto- 
nian which consists of the intrinsic kinetic energy and 
the Vucom interaction derived from the Argonne V18 
potential including the Coulomb potential, 

'H = Hint = T — Tcm + UUCOM = Tint + UuCOM , (5) 

in two-body approximation. The intrinsic kinetic energy 
operator reads. 

Tint = 7^ - T^cm = I ^ ^ (P, - )^ (6) 

Kj 

where we assume equal proton and neutron masses, m. 
It is the two-body Hamiltonian Hint of Eq. Q that has 
been used in Hartree-Fock (HF), perturbation-theory, 
and RPA calculations in Refs. |1(1 [llj and that will be 
employed in this work too. In practice, two-body ma- 
trix elements in a harmonic-oscillator basis are the in- 
put to such calculations. The harmonic-oscillator basis 
is typically restricted to 13 major shells, which warrants 
complete convergence of the HF results. Calculations 
with larger basis sizes are possible but increasingly time- 
consuming. The optimal value of the harmonic oscillator 
parameter is determined from an explicit energy mini- 
mization for different regions in the nuclide chart. 
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III. THE ERPA METHOD 

The method used in this work is a first-order RPA. 
The excited nuclear states are described as hnear combi- 
nations of particle-hole {ph) and hole-particle (hp) con- 
figurations and created by a set of vibration creation op- 
erators. The RPA ground state (reference state) is the 
vacuum of the corresponding annihilation operators and 
contains induced long-range correlations. As mentioned 
in Sec. in order to facilitate the actual solution of the 
RPA equations, one usually approximates the reference 
state by the uncorrelated HF ground state when calcu- 
lating the RPA matrix elements and the transition ma- 
trix elements. Here we use a renormalized formulation 
of the RPA that avoids this approximation and allows 
us not only to assess the effect of ground-state correla- 
tions on the excitation spectra, but also to evaluate some 
properties of the correlated RPA ground state. Following 
Ref. [Ill, we call this formulation Extended RPA (ERPA) 
to distinguish it from the one which is based on the un- 
correlated HF ground state. 

Next, we review the ERPA method based on Refs. 0, 
0, 01 and we present the formalism used to evaluate 
ground-state properties and transition strengths. Ex- 
pressions with explicit angular-momentum coupling are 
shown. Closed-shell nuclei are considered and spherical 
symmetry is assumed. The symbol p{h) will represent 
all the quantum numbers of a particle (hole) state except 
the magnetic quantum number TOp(m/i), i.e., the set of 
quantum numbers {n-p(/t)^p(?i)Jp(/i)ip(?i)} of the nodes, or- 
bital angular momentum, total angular momentum, and 
isospin. The letters a, f3, ... will be used to denote single 
particle states of either kind {p or h). As a single-particle 
basis, the one that diagonalizes the one-body density ma- 
trix (OBDM) is chosen — i.e., the natural-orbital basis. 
The ERPA formalism is based on this choice. 



A. The ERPA equations 

The creation operator Q\ of an excited state \v) 
with multipolarity J and parity tt is written as a linear 
combination of renormalized ph creation and annihilation 

JM 
ph 



operators A'^^f^, A:'^'^ 



ph 



i-^ph ■^ph 



Y;;;'\-iy+''A'^^''m,{7) 



where |0) is the ground state, for which 
Q.,j-A/|0) =0, 



(8) 



and the sum runs over ph pairs which can couple to to- 
tal spin and parity . The renormalized operators are 
written as linear combinations of ph operators, 

,.JM-\ _ nripjptp-Zhjhthr i iJM 



where the angular-momentum coupled operators are 
given by 

KahY""' ^ {jpmpjhmh\JM){-lY'+'^-al^^ah-m,. 



niprrih 



(10) 

The choice of the single particle basis and the require- 
ment that 



(11) 



following from the fermionic character of the operators, 
dictate that the matrix N be given by |3 



N 



npUh.n'pn'f^ ~ 



Snpn'Jn^n'^iPh - Pp) 
ph 



= ^npU'Jn^n'^D^^''^. (12) 



Here 



l4m„aam„|0) 



(13) 



denotes the occupation probability of the single-particle 
state a, which is independent of TOq.. It can be shown 
that, under the above circumstances, the amplitudes X 
and Y of Eq. |(7|) obey the orthonormalization condition 

E(<'''*^;r - Yph'^*Y;f) - '5..'. (14) 

ph 

The equations determining the amplitudes X and 
Y are obtained by using the equations-of-motion 
method . As usual, they are written as a matrix equa- 
tion in the ph space, 



A B 

-B* -A* 



X" 



X" 



V \ yv \ I 



(15) 



where E^, is the excitation energy of the state \v) and 
the quantum numbers J'^ are implicit. The elements of 
the matrices A and -B, which are independent of M, are 
given by 



php'h' 



and 



B 



p,p,,, = -mAii:',n,i-iy+"Air,f'] 



(16) 



(17) 



where Ti is the Hamiltonian of the system. The sym- 
metrized double commutators are defined as 

[A,B,C]^^{[A, [B,C]] + [[A,B],C]). 
Here we consider the two-body Hamiltonian, 

Ti. = -ffint = ^ X! Hai3a'i3'al^a^/jap>aa', (18) 

containing a two-nucleon interaction and the A- 
dependent intrinsic part of the kinetic energy, as intro- 
duced in Secini Eq. ©. 
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The remaining task is to evaluate the expectation val- 
ues of the double commutators in Eqs. l(T^ and ((T7|l . 
In general, elements of both the one-body and the two- 
body density matrices will appear in the commutator 
of the Hamiltonian iJint with a ph operator. The ele- 
ments of the two-body density matrix can be tackled us- 
ing a linearization procedure [Tol |2C| , by which they are 
contracted with respect to a reference state (the ground 
state) so that a quantity linear in the ph operators is ob- 
tained. Then, by substituting in Eqs. and ((T7|) . the 
matrix elements of A and B can be expressed in terms 
of the OBDM only. In the standard RPA, one uses the 
HF state (for which p/j = 1, = 0) as the reference 
state. In the present ERPA the correlated ground state 
0) is chosen as the reference state. By performing this 
linearization, one obtains 



A-' 



php'h' 



1 

2 ^ ph ^ p'h' 
X {Shh'£p'p ~ <5, 



d: 



-1/2^1/2 . 
ph p' h' ) 

pp' ^hh' ) 



-D 



1/2^1/2 

ph 



Dlil{ph-^;J\Hn,t\p'h' (19) 



[DphDp,h'{l + 5hh'){l 
x{{ph-'-J)ip'h' 



Spp')]'^' 



J-i)|i7int|0) (20) 



where \a(3; J) are the normalized and antisymmetrized, 
angular-momentum-coupled two-body states and 

=Y,Pp[{l + <5a/3)(l + 5^,[3)]'/^{a(3; J|i?int|a'/?; J) 

(21) 

are matrix elements of the single-particle Hamiltonian 
(see Sec. IIII B(l . The quantity Dpfi has been defined in 



Eq. (|12|l . For ph = 1, Pp = the usual RPA equations 
are recovered. 

According to Eqs. (HHJl and (j^OJl, the ERPA matrix ifT^ 
depends on the occupation probabilities pa- These, as we 
will see, depend on the forward and backward amplitudes 
X and Y of all the J'^ states of the nucleus, i.e., on the so- 
lutions of the ERPA equations. Therefore, the equations 
must be solved iteratively. Let us call this procedure 
the ERPA iteration. The occupation probabilities of the 
single-particle states can be expressed in terms of the am- 
plitudes X and Y and the matrix TV of Eq. (|12|l . using 
the number-operator method . Contrary to the quasi- 
boson approximation, this method does not suffer from 
double-counting the second-order contributions fill IT^ . 
For the OBDM in a generic basis one finds 0, Il5| 

(0|a],™pVmp,|0) = 

hirihi^i^' pihi 

X (0|4mpahmJi')(j^'|aL,^ap'm^,|0), (22) 

prripi/y' pihi 

X (0|at„^ah.,„ji/)(j/'|aj;„^^ap'„^,|0). (23) 

In the natural-orbital basis, and using angular momen- 
tum coupling, we have 



Pi hi 



uy'p 



Pihi 



(24) 



(25) 



where j = \J (2j -I- 1). These equations, which are exact 
up to terms 0(|y|*), are not linear in the OBDM ele- 
ments. Therefore, for given X and Y amplitudes (at a 
given step of the ERPA iteration), they have to be solved 
iteratively as well: one begins with an initial guess for 
the occupation probabilities on the r.h.s., calculates the 
OBDM on the l.h.s., diagonalizes it, thus obtaining a new 
set of occupation probabilities, and so on, until conver- 
gence is reached. Let us call this the OBDM iteration. 

This double iterative procedure of solving the ERPA 
equations can be outlined as follows: In order to begin the 
ERPA iteration, we need an initial choice for the ground 
state. We start with the HF ground state. We build 



the ERPA matrix for each J'^ that we wish to take into 
account by using Eqs. (|19|l and (|20|l and assuming the 
HF occupation probabilities (1 for holes and for parti- 
cles). Then we solve the eigenvalue problem of Eq. I|15|l . 
At this stage, we have simply solved conventional RPA 
equations. Then we calculate a "new" OBDM by using 
Eqs. H24|) and (|25|l and, as input, the X and Y ampli- 
tudes and the HF occupation probabilities. We proceed 
with a OBDM iteration and obtain new occupation prob- 
abilities, which hereafter replace the initial ones (the HF 
ones in this first step) and are used as input to build 
again the ERPA matrix for each J'^ . Here begins the 
second ERPA-iteration step. We proceed as before, until 
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convergence is reached. 



B. Ground-state properties in the ERPA method 

The ERPA method provides information on the corre- 
lated RPA ground state. As we have seen, the natural 
orbitals can be determined and their occupation proba- 
bilities can be calculated up to fourth order in the back- 
ward amplitude. As a convenient way to quantify the 
ground-state correlations we will use the mean square 
deviation per particle of the OBDM , defined by 



a-Tr[(e-^;(0))2]/A 



(26) 



This quantity characterizes the deviation of the OBDM 



Q of the correlated ground state from the OBDM g' 



(0) 

of the uncorrelated state described by a Slater determi- 
nant I '5*^")). It describes how well |\E'*^°)) approximates 
the correlated state. It is minimal when evaluated in 
the natural-orbital basis, with |\E''^°)) built of the A hole 
states l22l. Then 



1 

A 



. h 



E--2 2 
JpPp 



(27) 



expresses the depletion of the hole states and the occu- 
pation of the particle states as an average. In HF the- 
ory (Tmin obviously vanishes, as does the Fermi-sea de- 
pletion, while in realistic cases Omin is expected to take 
values of the order 10^^ [23, corresponding to a 
Fermi-sea depletion of at least 10%. It has long been 
observed [23| that values way below 10~^ indicate that 
short- and medium-range correlations have not been (ad- 
equately) taken into account by the model in question. 
This is typically the case, e.g., in RPA calculations. We 
will return to this point in Sec. IIVI 

Removal and pick-up centroid energies, e^"-* and e^"*"-*, 
respectively, can also be defined 18, 20, 2^ and evaluated 
within the ERPA. The removal energy e"^^ correspond- 
ing to the hole state h is the centroid of the distribution 

= {0\ai^J{E-H,^,)ah^jO), (28) 

where |/), Ef are the eigenstates and energies of the 
{A — l)-system, i/i„t|/) = Ef\f). One finds, using again 
a linearization procedure, 

4"^ ^ -J2pM + ^c>h){ah;J\Hint\ah;J) = -Sh- (29) 

a, J 

Similarly one finds the centroid of the energy correspond- 
ing to picking up a particle at the state p, 

e^+^ = ^Pq(1 + 5ap){ap; J\Hint\a.p; J) = Sp. (30) 

a, J 



Defined in the natural-orbital basis. Eh and Ep no longer 
correspond to the eigenvalues of the single-particle (s.p.) 
Hamiltonian — in contrast to the HF case. A s.p. Hamil- 
tonian can be defined, in principle, as the one-body 
term of the nuclear Hamiltonian operator, Eq. H18I) . when 
written in a normal-ordered form 20, 25] . It differs from 
the HF s.p. Hamiltonian only by a similarity transfor- 
mation and therefore posesses the same eigenvalues. If, 
however, one linearizes the two-body part of the normal- 
ordered total Hamiltonian iJint and adds the thus ob- 
tained one-body operator to the genuine one-body part 
i/jjj^ , one obtains a s.p. Hamiltonian Hl^^' whose ma- 
trix elements in the natural-orbital basis are given by 
Eq. (|21|) . Its diagonal matrix elements coincide with the 
Eh, Ep energies defined in Eqs. (j^ . It is possible 

to diagonalize the s.p. Hamiltonian in order to obtain its 
eigenvalues, Cq. 

All s.p. energies defined above have been formally 
evaluated in terms of matrix elements of the two-body 
Hamiltonian i/int. They should be further corrected for 
the center-of-mass energy, as was done in Ref. for 
the HF s.p. energies. For the purposes of this work, 
however, it will suffice to compare our uncorrected re- 
sults with the uncorrected HF energies, in order to show 
the relative effect of RPA correlations. 

One should keep in mind, that only for s.p. levels 
around the Fermi energy is it sensible to discuss discrete 
s.p. energies. For low- lying holes or high- lying particles 
the s.p. strength is strongly fragmented in realistic cases 
— mainly due to short-range correlations j2^. 



C. Transition strength 

We consider the response of the nucleus to an external 
field described by the s.p. multipole operator 

Oj^M= ^ {ama\Oj^M\f3mf3)a''i^^^aam^. (31) 
The corresponding strength distribution is given by 



i?(i?)^^|(Ho, 



H{E^E,), (32) 



where the sum runs over states with quantum numbers 
. The transition matrix element between the ground 
state and the excited state [see Eq. Q] is given by 

{v\0.j.m\Q) = 

(ama|Oj^M|/?m^) 



X (0|Q,,j.A/a]j,„^a„„jO). (33) 
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After we substitute for the operator Q using Eqs. |7J, ((HI, 
and H12|l we obtain 

{v\Oj.m\Q) = 

ph 

>^Dli^{p\\Oj\\h). (34) 
The s.p. transition ampUtudes are renormahzed by a 

1 /2 

factor due to the ground-state correlations. 

In the UCOM framework the operators of all observ- 
ables must be transformed in a consistent way. In par- 
ticular, the same unitary transformation as used for the 
nuclear Hamiltonian has to be applied to the transition 
operators as well. Thus, correlated transition operators 
C'^Oj^C should be used in Eq. However, it has 

been shown that the effect is negligible, at least for 
the monopole and quadrupole operators. 

Having solved the ERPA equations, one can use 
Eqs. H32|l and H34|) to calculate the response of nuclei to 
external fields of interest. In this work we examine the 
response to isoscalar (IS) and isovector (IV) operators 
with natural parity, tt — (—1)'^, of the usual form 

A 

1=1 

A 

Oy^M = T.^P^^'YjMin) (36) 

i=l 

for J ^ 1, with = J for J > 2 and = 2 for J = 0. 
For the dipole operators we will use the corresponding ef- 
fective forms which remove the mixing with the spurious 
center-of-mass motion, 

0?-M = E(^' - ^{rfn)YMr^) (37) 

^=l 

Z ^ N ^ 

O'^-M = ^Y.'-^Y,M{f{)-^Y.'^,Y^^t{h). (38) 

1=1 1=1 

Thus, the spurious state is automatically eliminated from 
the strength distributions. 

Allowing for non-zero depletion (occupation) of the 
hole (particle) states in the ground state means that hole- 
hole and particle-particle transitions should be allowed, 
besides the ph ones. An extension of the presented ERPA 
scheme has been introduced j23| to take these additional 
configurations into account. It was shown that, due to 
their omission, the energy-weighted sum rules are not 
preserved by the ERPA. For closed-shell nuclei, however, 
like those considered here, the deviation should be small. 

IV. RESULTS 

We have applied the ERPA method to examine the re- 
sponse of closed-shell nuclei to isoscalar (IS) and isovec- 
tor (IV), natural-parity external fields. We will present 



results obtained using the correlated Argonne V18 inter- 
action with 1^^'°'' = 0.09 fm^. In most cases we have used 
a harmonic-oscillator basis of 13 shells, with an oscillator 
parameter ano = 1-8 fm. We will call this set of param- 
eters "standard", for brevity. The term will mean also 
that both natural- and unnatural-parity modes (but not 
charge-exchange ones) of J < Jmax = 4 have been taken 
into account in the evaluation of g.s. correlations. In 
certain instances we have changed the size of the basis, 
to test the convergence of our results, or used a different 
"/max value, or only natural-parity states. 

The main question is how our RPA results, such 
as those presented in Ref. will be affected by 

the ground-state correlations taken into account by the 
present scheme. It is instructive to first discuss the 
ground-state properties calculated within ERPA. Indeed, 
the s.p. occupation probabilities enter directly both the 
ERPA equations, Eqs. (UU, (jTHI), HOI), and the calcula- 
tion of the transition strength, Eq. H34|) . They renormal- 
ize the matrix elements of the interaction in the former 
case and of the ph transition operator in the latter, via 
the quantities Dph- For given s.p. states, a more depleted 
Fermi sea would imply a weaker residual interaction and 
reduced ph transition strength. Of course, the s.p. eigen- 
states are changed by the self-consistent solution of the 
ERPA equations and affect the final result as well. 



A. Ground-state properties 

For the nuclei ^^O, ^Ca, ^"Zr, and ^""Sn and for the 
standard parameter set defined above we have evaluated 
the HE eigenenergies, the eigenenergies Cq. of the new 
s.p. Hamiltonian and the centroid energies defined in 
Eqs. and l|Sn|) . None of these sets of s.p. energies 
are corrected for the center-of-mass energy. 

In Fig. n we show our results for the ^''Zr nucleus. 
The thick grey line separates the particle states from the 
hole states. One can observe that the spectrum of the 
new eigenenergies Cq, is more dense than the HF one. 
The centroid energies £p and et are even closer to each 
other, on average, while the Fermi gap in this case ap- 
pears larger. Indeed, the lowest pickup centroid energy 
corresponds to the strength distribution to one or more 
particle eigenstates — since there is no mixing between 
particle and hole configurations. Therefore, its value can- 
not be lower than the lowest particle state's Cq. Similar 
arguments explain why Bq > for the highest hole level 
and Cq < Ea for the most bound one. Let us note also 
that, in general, the ordering of the levels is different 
than for the Cq and HF ones. In total, the ph energies 
Cp — eh are smaller than their HF counterparts which 
enter the RPA equations. This means that the "unper- 
turbed spectrum" corresponding to the ERPA equations 
(the one obtained by setting the B matrix and the last 
term in Eq. H19|l equal to zero, and assuming that the 
prefactor of the first term in Eq. (|19|l is roughly equal to 
one) will be denser than the unperturbed HF spectrum. 
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FIG. 1: (Color online) The Hartree-Fock eigenenergies, the 
eigenenergies Ca of the new s.p. Hamiltonian and the pick-up 
and removal centroid energies ±ea, Eqs. for the 

nucleus ^"Zr, evaluated using the standard parameter set (see 
text). The energies are not corrected for the center-of-mass 
energy. A thick grey line separates the particle states from 
the hole states. 
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The above remarks hold for all examined nuclei. 

In Fig. 21 the deviation Apa = Pa — Pa'' of the occu- 
pation probabilities pa from their HF values pi"^ (1 for 
holes, for particles) is shown for the proton s.p. levels 
of the nuclei ^'^Ca and ^"Zr. The full circles correspond 
to the standard parameter set. In the case of the open 
circles, only natural-parity excitations have been taken 
into account. More than half the total depletion seems 
to be accounted for by these excitations. 

In Table we list the values of the mean square de- 
viation per particle cr, Eq. 1)2 7|) . that we find for the 
four nuclei ^^O, "OCa, ^''Zr, ^''OSn when all J± excita- 



TABLE I: The values of the mean square deviation per 
particle a, Eq. 1271 . found when all J* excitations (itnup 
- standard parameter set) and only natural-parity (ctnp) ones 
are considered. In square brackets, for '"Zr, results obtained 
within a larger s.p. basis (15 shells) are shown. 





16 Q 


«'Ca 


*Zr 


lOOSn 


fJNUP X 10^ 


3.43 


2.23 


1.11[1.17] 


1.37 


CTNP X 10^ 


1.85 


1.05 


0.46[0.48] 


0.57 



FIG. 2: The deviation Ap^ = pa~pa of the occupation prob- 
abilities Pa from their HF values pi"' is shown for the proton 
single-particle levels of the nuclei *°Ca and ®°Zr. Filled cir- 
cles: standard parameter set; open circles: only natural-parity 
states taken into account. The energies are not corrected for 
the center-of-mass motion. 



tions (ctnup - standard parameter set) or only natural- 
parity ones (cnp) are considered. (Note that no charge- 
exchange excitations are considered.) a characterizes the 
depletion of the Fermi sea due to correlations. In square 
brackets, for ^"Zr, results obtained within a larger s.p. 
basis (15 shells) are shown, to verify that our conclusions 
remain valid. The ratio (Tnp : ctnup is roughly 50%. 
We notice that the values of a are of the order 10~^, 
reflecting the effect of long-range correlations. As men- 
tioned in Sec. IIII Bl such cr-values are typical of models 
which do not take into account short-range correlations. 
In principle, these are dealt with by the UCOM, but a 
full treatment to this end has not been pursued here. 
Indeed, when calculating the OBDM and the occupation 
probabilities pa , we have used the uncorrelated operators 
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ajjOo,/. For a consistent UCOM treatment, however, cor- 
related operators C^al^aa'C should be considered. Such 
calculations are not trivial and lie beyond the scope of 
the present work. As a relevant example where the po- 
tential of the UCOM was fully exploited in describing 
short-range correlations in nuclei, we mention the calcu- 
lations of Ref. 8] , where the high- momentum tail of the 
s.p. momentum distribution is reproduced. 

For "^^Ca and ^°Zr we can compare our present results 
for the depletion of the Fermi sea to the ones obtained 
within second-order perturbation theory (PT) and re- 
ported in Ref. In general the occupation proba- 
bilities of hole states are larger in the case of ERPA, 
especially for ^'^Zr. We note that within PT more con- 
figurations are taken into account in the evaluation of 
ground-state correlations, namely all the values of 
available within the basis, including charge-exchange ex- 
citations. More importantly, the use of non-renormalized 
PT is known to overestimate the depletion of the Fermi 
sea [2^. This can explain also why for ^°Zr the parti- 
cle occupation probabilities obtained within ERPA are 
much smaller. For "'"Ca they are larger only close to the 
Fermi energy. For both nuclei the ERPA values of the 
particle occupation probabilities drop much faster to zero 
as the energy increases, compared to PT. The use of the 
natural-orbital basis, rather than the HF one, is respon- 
sible for this behaviour; it is not a matter of RPA vs. PT 
— cf. the particle occupation probabilities reported in 
Ref. m within (S)RPA in the HF basis. 



B. Collective excitations 

In Fig. Owe show the centroid energies of the IS GMR, 
IV GDR, and IS GQR obtained for the nuclei i^Q, 40Ca, 
^•^Zr, ^°°Sn, and ^°^Ph within the RPA (solid bars) and 
ERPA (dotted bars), compared with some experimental 
data (points). In Fig. 21 we show the energy of the low- 
lying IS octupole (ISO) state. As examples, we show 
in Figs. and the corresponding IS monopole, IV 
dipole, and IS quadrupole strength distributions for -'^^O 
and ^"Zr. Also shown in these figures are the IS dipole 
strength distributions. In all cases shown in Figs. 13 -0 
the standard parameter set was used, except for ^°^Pb, 
where a basis of 15 shells was employed to ensure bet- 
ter convergence in this heavy nucleus and Jmax = 3 was 
used to avoid an extremely time-consuming calculation. 
The strength distributions up to 50 MeV were taken into 
account when calculating the centroid energies, which 
are defined as the first energy moment of the distribu- 
tion, mi, divided by the total strength, toq- The ex- 
perimental centroids of the IS GMR and the IS GQR 
were taken from Refs. JOJ (^^O), C'°Ca), 32] (^"Zr), 
and (^°®Pb). Photoabsorption cross sections were 
found in Refs. Ip (^^O), [H ('^''Ca), |36i (^^Zr), and 
[33, nil (^°*Pb) and the centroids toi/toq of the corre- 
sponding IV GDR strength distributions were evaluated 
from those. Experimental energies of the collective oc- 



tupole state were adopted from Ref. |39|. 

From Fig. 13 it is evident that the energy of the GDR 
drops when ground-state correlations are explicitly con- 
sidered within the ERPA, but not enough to reach the 
experimental data. In general, a decrease of no more 
than 1 MeV was achieved for the heavier nuclei. The 
same holds for the IV 0"*" and 2+ centroids (not shown). 
The IS giant resonances, namely the IS GMR, GDR, and 
GQR, were less affected. In most cases, their energies 
were higher when evaluated within ERPA than within 
RPA. Apparently, in the IS case the weakening of the (at- 
tractive) residual interaction due to the depleted Fermi 
sea compensates for the compression of the s.p. spectra. 

The energy of the collective IS octupole state. Fig. 01 
increases by up to 1 MeV and the disagreement with 
the experimental data is worsened. We found that its 
strength within ERPA decreases with respect to RPA 
by 10-20%. Still, the ERPA values of the strength are 
larger than the experimental ones by a factor 1.6-2.2 for 
the nuclei considered. The energy of the low- lying IS 
quadrupole state for the nuclei ^°Zr, ^^'^Sn, and ^ospj^ 
(not shown) increases due to the correlations by no more 
than 0.3 MeV and its strength decreases by 10-20%. 

We found that, for all nuclei and excitation fields ex- 
amined, the total strength ttiq and total energy-weighted 
strength mi decrease within ERPA. The relative change 
remains below 10% in almost all cases, being largest for 
the lightest nuclei. Within ERPA the spurious dipole 
state is found at approximately 1 MeV, i.e., somewhat 
further from zero than within RPA. 

In Fig. [71 we show, for ""^Ca and ^''Zr, the cumula- 
tive effect on the IS GMR and IV GDR centroids of al- 
lowing zero-point fluctuations of ever larger multipolar- 
ity Jmax- The crosses mark results obtained using only 
natural-parity states. The squares mark results obtained 
using both natural- and unnatural-parity states. Lines 
are drawn to guide the eye. The horizontal lines show the 
RPA results for comparison. In calculating the centroids, 
the corresponding strength distributions up to 50 MeV 
were taken into account. For ^"Zr and Jmax = 4 an ex- 
tra set of data is shown in each panel, detached from the 
lines. These were obtained using a larger s.p. basis of 15 
shells and verify that our calculations have converged to 
a rather satisfactory degree with respect to the basis. 

In all cases a steeper change is observed between 
Jmax — 2 and 3, indicating the important effect of the 
low-lying collective 3^ state. We see also that the conver- 
gence of our results as Jmax increases is rather slow. As 
RPA calculations of ground-state energies suggest , in 
order to achieve good convergence when describing RPA 
correlations one may need to take into account Jmax val- 
ues up to approximately 10. For ERPA, such large val- 
ues are beyond our computational capabilities at present. 
Since, however, no strongly collective states with multi- 
polarities larger than 4 are expected, and based also on 
the results of Ref. |3 , it is reasonable to speculate that 
the difference between the RPA and ERPA results would 
be no more than doubled if higher values of Jmax were 
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FIG. 3; (Color online) Centroid energies of the IS GMR, IV GDR, and IS GQR obtained for the nuclei ^^O, ''°Ca, *Zr, ^°°Sn, 
and ■^"^Pb within the RPA (red solid bars) and ERPA (blue dotted bars), compared with experimental data (black points — 
for references see text). The standard parameter set was used, except for ^'^*Pb, where Jmax = 3 and a s.p. basis of 15 shells 
were employed. In calculating the centroids, the strength up to 50 MeV was taken into account. 
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FIG. 4: (Color online) As in Fig. |21 energy of the low-lying 
IS octupole (ISO) state. 



considered. Of course, no charge-exchange excitations 
have been taken into account here. These can introduce 
additional correlations of non-negligible ampHtude, albeit 
smaller than the non-exchange ones ,12] . 



V. SUMMARY AND PERSPECTIVES 

A correlated realistic interaction, derived from the Ar- 
gonne V18 nucleon-nucleon interaction in the UCOM 
framework, was employed in calculations within the so- 
called Extended RPA (ERPA). The aim was to inves- 
tigate to which extent such an interaction can describe 
nuclear collective motion in the framework of first-order 
RPA. Renormalized formulations of the RPA, such as the 



ERPA, consider explicitly the depletion of the Fermi sea 
in the ground state due to RPA correlations. In a pre- 
vious work we used the same interaction within a RPA 
model [l]| . where excitations are built on top of the un- 
correlated HP ground state. The ERPA employed here 
allowed us to examine the effect of ground-state correla- 
tions on the excitation spectra. It was found that the ef- 
fect on the properties of giant resonances is rather small. 
Compared to the standard RPA model, their centroid en- 
ergies decrease by up to 1 MeV in the isovector channel. 
The isoscalar response is less affected in general. Ground- 
state properties obtained within the ERPA were com- 
pared with corresponding HF and perturbation-theory 
results and discussed as well. 

Our results with the correlated Argonne V18 show that 
the RPA based on the uncorrelated HF ground state is 
a fairly good approximation to the actual RPA prob- 
lem in the case of giant resonances, although with this 
interaction the HF solution is a bad approximation to 
the RPA ground state when one considers the nuclear 
binding energies. It should be noted that a theoreti- 
cal error of 1 MeV in the centroid energies (or more, 
if charge- exchange and highcr-mutipole correlations are 
included) can be important in certain applications. That 
said, the inability of our (E)RPA calculations to repro- 
duce the experimental data should originate elsewhere. 
Up to now we have assumed that residual three-body 
forces can be neglected, based upon the fact that they 
contribute only marginally to the ground-state energy as 
calculated within many-body perturbation theory |l7j . 
This is not necessarily a valid assumption. We are cur- 
rently constructing a simple phenomenological zero-range 
three-body force, to be used along with the correlated 
two-nucleon interaction in future calculations. Prelimi- 
nary results show that by using such a three-body force 
it is possible to improve on the description of observables 
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FIG. 5: (Color online) IS monopole, IS quadrupole, and IS and IV dipole strength distributions of ^^O computed within RPA 
(red solid lines) and ERPA (blue dotted lines). The discrete strength distributions have been folded with a Lorentzian of width 
r = 2 MeV. The standard parameter set was used. 




FIG. 6; (Color online) As in Fig.E) for ""Zr. 
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FIG. 7: (Color online) The cumulative effect of allowing zero-point fluctuations of larger multipolarity Jmax on the IS GMR 
and IV GDR centroids of '"'Ca and '"Zr. Blue crosses: using only natural-parity states. Red squares: both natural- and 
unnatural-parity states were used. A s.p. basis of 13 shells was used. Lines are drawn to guide the eye. The horizontal lines 
show the RPA results for comparison. For ®°Zr and Jmax = 4 an extra set of data is shown in each panel, detached from the 
lines. These were obtained using a larger s.p. basis of 15 shells. 



such as nuclear radii and resonance energies while retain- 
ing the good reproduction of the experimental binding 
energies. 

Another important issue is that, within the RPA, 
one neglects the coupling to higher-order configurations 
{2p2h and beyond). Collective excitations can be signif- 
icantly affected by the inclusion of 2p2h configurations, 
as preliminary SRPA results show. Therefore, SRPA cal- 
culations will also be the subject of future work. In prin- 
ciple, it is possible to combine the SRPA with a corre- 



lated ground state |40|,|41j for a most complete theoretical 
treatment of nuclear excitations. 
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